Kramers equation for a charged Brownian particle: The exact solution 
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We report the exact fundamental solution for Kramers 
equation associated to a brownian gas of charged particles, 
under the influence of homogeneous (spatially uniform) oth- 
erwise arbitrary, external mechanical, electrical and mag- 
netic fields. Some applications are presented, namely the 
hydrothermodynamical picture for Brownian motion in the 
long time regime. 



I. INTRODUCTION 

In Chandrasekhar's 1943 celebrated paper 1 , Kramers 
equation 2 was solved for the free Brownian particle and 
some general lines were drawn towards solving this prob- 
lem in a field of force. Only recently some progress 
was reported considering Kramers equation for a charged 
Brownian particle in a field of force: Czopnik & Gar- 
baczewsky (CG) 3 solved Kramers planar equation in 
a magnetic field, essentially transforming the magnetic 
field contribution into a tensorial Stokes-like dissipative 
force. Later, Ferrari 4 via transformed phase space vari- 
ables, mapped Kramers equation for a charged Brownian 
particle in an electric field into the free Brownian particle 
case. By combining both CG's 'rotated' Stokes force and 
Ferrari's gauge, in section II we report the exact funda- 
mental solution of Kramers equation for a charged Brow- 
nian particle in an uniform, otherwise arbitrary field of 
forces. In section III we present some applications, con- 
cerning the hydrothermodynamical picture of Brownian 
motion, the validity of the local equilibrium approxima- 
tion and the 'linear' regime (see for example 5 ' 6 ). Com- 
parison is made with some results obtained via a per- 
turbative recursive scheme 7 ' 8 . Finally in section IV we 
present some concluding remarks and outline some work 
in progress. 



II. FUNDAMENTAL SOLUTION 

We study a Brownian gas composed of charged parti- 
cles (mass m, charge e) under the influence of external 



fields: mechanic (mec), electric (E) and magnetic (B) 
fields, uniform in space and in general time dependent. 
Our starting point is Kramers equation for the density 
probability distribution P(x, v,i) in phase space (posi- 
tion r, velocity v) at time t, in contact with a reservoir 
at temperature Tr and under the force fields 



F(v,t) - F mec + eE + -v x B 
The associated Kramers equation 2 reads 



(1) 



dP 
~dt 



dP 
<9x 



F dP 

m dv 



d F d T R d 2 P 
ov m m av 2 



with Boltzmann's constant jfcg = 1, A = t _1 is the 
friction coefficient (inverse of the collision time) and 
Fd = — Amv is the dissipative Stokes-like force. As in 
previous work 7 ' 8 , where a perturbative recursive scheme 
was presented for a more general case, we define vt, a 
thermal velocity given by mv T = Tr and dimensionless 
variables by scaling space, velocity and time, respectively 
with I = tvt, vt and r. Also we define the conservative 
acceleration a (and the associated potential cf>) and the 
cyclotronic frequency vector $7 respectively as (hereafter 
all quantities are dimensionless unless stated otherwise) 



(Fr, 



mvT 
SJ = B =uj u) 



eE) = --^ 
ax 



(3) 
(4) 



Notice that in terms of the usual dimensional cyclotronic 
frequency 8 , we have u> — uj c t. Concerning notation, we 
chose the very convenient bra-ket convention. Denote 
any vector V by V =V X \x)+V y \y) + V z \z) and its adjoint 
by V T = V x (x\ + V y (y\ + V z (z\ (all quantities are real). 
We also define some useful dyadics, namely: ei = \z) (z\, 
e 2 = \x) (x\ + \y) (y\, e 3 = \x) (y \ - \y) (x\ and the unit 
dyadic e = ei + e2. Furthermore, we define the z-axis as 

the magnetic field direction (u>=z), and Kramers equa- 
tion (2) is cast in a compact form as 

8P dP dP d A „ d 2 P 
at ax ov av ov z 

where the magnetic contribution is included as a tensorial 
Stokes-like dissipative term (see 3 ' 8 ) 
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Av = v + fi x v =(e — u;e 3 )v 

We also define 

M = A _1 = ei + a (e 2 + cue 3 ) 

with a -1 = 1 + uj 2 . In dimensional form this tensorial 
collision time constant has the familiar form 8 



1 U) r T 



M 



-LO c T 



1 + (c c r)2 I Q Q 1 + (wcr) 



(6) 



Since the planar (a;, y) dynamics is decoupled from the 
z axis dynamics, the case a = can be trivially solved 
generalizing the planar results of (CG) 3 . Then, trans- 
forming to new variables R and V, via Ferrari's gauge 4 
we map our problem to the free Brownian particle solved 
in 1943 by Chandrasckhar 1 . Here we particularize to 
time independent external fields case, yielding Ferrari's 
transformed variables 4 

R = x (5x(a, M,t) 
V = v - tfv(a, M,t) 

where 

5\ = M(l - 0)a + 0v o 

5x = Mat - M 2 (1 - 0) a + M (1 - 0) v + x 

= exp(— t) (e 1 + e 2 cos ut + e 3 sin ut) 

and with Kramers equation (5) mapped to the trivially 
solved 3 ' 4,1 equation 

dP flP d ATrn d 2 P 
h V = AVP H 

dt + or dv ^ dv 2 

The fundamental solution G(x, t,i|x , v ), namely with 
free boundary conditions and the point-like initial condi- 
tion 

G(x, v,i = 0|x , v ) = S(x - x )5(v - v ) 
is given by 

1 1 



G(x,t,i|x ,v ) 



== cxp — -<l? 



(7) 



where 



with 



* d = V f A,,V + R f A r R 
* sn = VtA ro R + RtA m V 
* a = 2Qt (R x V) 



* = ( Vt Rt ) A 



R 



where 



with 



and where 



A„ -A m + |Q|e 3 

-A m -|Q|e 3 A,. 



A fl " i a a 

A Q = ^e 2 + — ei 



fc A 

A = a„a r — a 2 ^ — A: 2 



a r = 1 — 6„ 



a m = a (l - 26 e 6 c + 6 2 ) 



fc = a (2b e b s — uoa r ) 



a v = a (a r + 2t — Aa (1 + 6 e (w6 s — 6 C ))) 



(8) 



b e = exp(-i), 



cos wi, 



sinwt 



The superscript * in equations (7, 8) denotes the corre- 
sponding quantity evaluated at null magnetic field (u = 
0). 

The general solution satisfying the initial condition 
P(x, v,t = 0) = Po(x, v ) is given by 

P(x, v,i)= y dx (ivoG(x, v,i|x , v )P (x , v ) (9) 

We normalize the probability P to TV, the total number 
of particles in the gas and define the particle density 



i(x,t) = J dvP(x,v,t) 



requiring 



or in a compact form 



N = J dxdvP (x,v) = J dxdvP(x,v,t) 

A mean density no is defined in the thermodynamic limit 
as N — uqV where the volume is defined by V — J dx. 
Thermal equilibrium conditions (TEC) are reached in 
the asymptotic regime t — > oo under null external fields, 
where we retrieve the Maxwellian distribution 



Ptec(v)= (2tt) 3/2 n exp--v 2 
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III. HYDROTHERMODYNAMICS OF 
BROWNIAN MOTION 

We start by defining some quantities of interest 7 ' 8 : par- 
ticle flow density J , the associated hydrodynamic veloc- 
ity u, kinetic energy density e, local gas temperature 9 
(in Tr units), local pressure p, entropy density s, Gibbs 
energy density g and the total chemical potential ^, re- 
spectively defined by 



/ 



J(x,t) = / rfvvP(x,v,t) = n(x,t)u(x,t) (10) 



■(x,i) = |n(x,t)0(x,i) = iy dvv 2 P(x,v,t) 



p(x,t) = n(x,t)9(x,t) 



s(x,t) = — J dvP(x, v, i) In xP(x, v, t) 



g(x, t) = e(x, t) - 0(x, t)s(x, i) + p(x, i) 



n(x,t) 



The additive constant 



g(x,t) 
n(x, i) 



lnx = -1 + 3 In 



tTr 



is chosen such that under TEC we retrieve the usual 
thermodynamical entropy density 7,8 . This equilibrium 
entropy and the associated chemical potential are respec- 
tively given by (in dimensional form) 



s eq (n ,T B ) = n - + In 



n Q (T R ) 



n 



Veq {n ,T R ) = -TrIii 



n Q (T R ) 



n 



with 



n Q (Tr) 



(2-KmT R \ 2 



\ h 2 J 

Furthermore we define an entropy flux density as 
J s = — J rfvvP(x, v, t) In xP(x, v, t) 

Balance equations are computed as in 7 ' 8 , yielding 
the continuity (Smoluchowsky) and the entropy balance 
equations, given respectively by 



and 



dn „, „ 



ds(x,t) dJ s (x,t) 

H ~ = °\ X >*J 



dt dx 

where the entropy production density a(x,t) is given by 



a(x,t) = J dvP(x,v,t) ^ 



<91nP(x,v,i) 
<9v 



3n(x,i) 



In this brief report we compute the non equilibrium 
long time regime (LTR, the limit t » 1 and with non 
zero external fields). In this limit we have: 

Sv = -M-^- = V F 
ox 

Sx = V F t + x* x* = xo - M(V F - v ) 

Thus, in the LTR Ferrari's 4 transformation is simply rep- 
resented by the constant velocity shift Vp or equiva- 
lently, given by the solution of the (dimensional) equation 

Vp 1 / e \ 

—=—[ F mec + eE + - V F x B 
t m \ c / 

We highlight some results in the LTR, corroborat- 
ing some previously obtained results, via a recursive 
method 7 ' 8 . We may cast 

P(x, v, t) = n(x, t)W(x, v, t) 

where to lowest order in t -1 we have 



Ant 



At 



n(x, t) 



1 \ 3 1 1 
— -exp--r(R) 
Ant J a F 2 ' 



with 



r(R) = ^(R 2 + a (fiR) 2 ) -V 2 D 



V D = -(R|Rxfl) 

Then it is also satisfied the relation 

, , ( dn dd> 
3 ^ = - M [8x +n £ 



(11) 



We define a magneto covariant derivative as in 8 for any 
given function /(x, t) as (notice the potential <fi is ex- 
pressed in Tr units) 
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D x f = Mexp(-^)— /(x,i)exp(0) 

and from (10) and (11) we retrieve Smoluchowsky equa- 
tion with a magnetic field 8 



dn 



d_ 

dx 



(D x n) 



(12) 



Also, the local temperature T(x, i), the entropy density 
and the total chemical potential (in dimensional form) 
are given by the following expressions 



T(x,i) 



= 6{x,t) = 1 



1 

2b 



1 



u 2 (x,t) 



(13) 



IV. CONCLUDING REMARKS 

Here we have presented the fundamental exact solu- 
tion for Kramers equation in a field of uniform forces, 
hitherto unknown, and applied the results to the long 
time regime. Work in progress seeks solutions for general 
initial conditions and other than free boundary condi- 
tions (for example membranes 10 ) and the associated hy- 
drothermodynamical picture, to be inserted into a more 
general framework 5 ' 6 . Also, we will address the ques- 
tion of a nonuniform reservoir temperature Tr(x) 8 ' 11 and 
the inclusion of chemical reactions 8 ' 12 to our Brownian 
scheme. 
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s(x, t) = s eq (n(x,t), T(x,i)) + Ss(x,t) 



(14) 



/x(x, b) = fi eq (n(x,i), T(x,t)) + 5fi(x, b) 



(15) 



where the entropy and chemical potential shifts, are re- 
spectively given by (to lowest order in b~ x and u 2 ~ M 2 ) 



Ss(x,t) = — -n(x,b)u 2 (x,b) 



5/i(x,t) - -T(x,i)u 2 (x,t) + 0(x,t) 

Furthermore, the hydrodynamic velocity u may be cast 
in several equivalent forms, namely 



u(x,i) — D x h\n — M- 



d ( fi(x,b) 



<9x \6(x,t) 



= V D +V F 



Lat us briefly comment on the last equation. The entropy 
density (equation 14) with u(x,t) = Z^lnn is formally 
reminiscent of a Ginzburg-Landau expansion as noted 
in 8 . Also, as stated by Landauer 9 the gradient of a bona 
fide nonequilbrium chemical must be proportional to the 
particle flux J = mi. 

It was found 8 that a suitable expansion parameter is 
the collision time r (or its tensorial partner M). In 
the LTR and confirming our conclusions in 8 , for the 
Brownian motion of a charged particle, the local equi- 
librium hypothesis is satisfied only bo firsb order (linear) 
in t(M) and where, among other things, Onsager rela- 
tions are satisfied. In this linear approximation we have 
T(x,b) w T R , Ss(x,b) w 0,6/j,(x,t) w <£(x,i). Furthermore 
up to second order in r 



If * l9< t> 



yielding a nonequilibrium temperature independent of 
the magnetic field, and a positive definite entropy pro- 
duction that vanishes as 
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